Long-lived dipolar molecules and Feshbach molecules in a 3D optical lattice 
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We have reahzed long-hved ground-state polar molecules in a 3D optical lattice, with a lifetime 
of up to 25 s, which is limited only by off-resonant scattering of the trapping light. Starting from 
a 2D optical lattice, we observe that the lifetime increases dramatically as a small lattice potential 
is added along the tube-shaped lattice traps. The 3D optical lattice also dramatically increases the 
lifetime for weakly bound Feshbach molecules. For a pure gas of Feshbach molecules, we observe a 
lifetime of >20 s in a 3D optical lattice; this represents a 100-fold improvement over previous results. 
This lifetime is also limited by off-resonant scattering, the rate of which is related to the size of the 
Feshbach molecule. Individually trapped Feshbach molecules in the 3D lattice can be converted to 
pairs of K and Rb atoms and back with nearly 100% efficiency. 

PACS numbers: 03.75.-b, 37.10.Pq, 67.85.-d, 33.20.-t 



Controllable long-range and anisotropic dipole-dipole 
interactions can enable novel applications of quantum 
gases in investigating strongly correlated many-body sys- 
tems [1-6]. Recent experiments have realized an ultra- 
cold gas of polar molecules in the ro-vibrational ground 
state [7] with high-resolution, single-state control at the 
level of hyperfine structure [8]. However, an obstacle to 
creating long-lived quantum gases of polar molecules was 
encountered with the observation of bimolecular chemi- 
cal reactions in the quantum regime [9]. Even with the 
demonstrated strong suppression of the reaction rate for 
spin-polarized fermionic KRb molecules, the lifetime of 
a 300 nK sample with a peak density of 10^^/cm^ was 
limited to ^1 s. Furthermore, when an external elec- 
tric field is applied to polarize the molecules in the lab 
frame, the attractive part of the dipole-dipole interac- 
tion dramatically increases the ultracold chemical reac- 
tion rate, reducing the lifetime of the dipolar gas to a few 
ms when the lab-frame molecular dipole moment reaches 
0.2 Debye [10]. A promising recent development was the 
demonstration that confining fermionic polar molecules 
in a ID optical lattice suppresses the rate of chemical re- 
actions even in the presence of dipolar interactions. Here, 
the spatial anisotropy of the dipolar interaction was ex- 
ploited by confining a gas of oriented KRb molecules in a 
two-dimensional geometry to suppress the attractive part 
of the dipolar interaction and thus achieve control of the 
stereodynamics of the bimolecular reactions [11]. In this 
regime, the lifetime of a trapped gas of polar molecules 
with a lab- frame dipole moment of ~0.2 Debye, a tem- 
perature of 800 nK, and a number density of 10^ cm~^, 
was ^1 s. 

In this letter, we study KRb molecules confined in 2D 
and 3D optical lattice traps, where we explore the effects 
of the lattice confinement on the lifetime of the ultracold 
gas. We note that a lifetime of 8 s has been achieved 
for homonuclear Cs2 molecules in a 3D lattice [12]. In 
our work, we find that long lifetimes are achieved for the 
molecules in a strong 3D lattice trap, even when there is a 
significant dipole moment in the lab frame. In addition. 



we observe that adding a weak axial corrugation to a 
2D lattice can result in long lifetimes for the trapped 
molecules. 

The experiments start with an ultracold mixture of 
2.9 X 10^ ^°K atoms and 2.3 x 10^ ^^Rb atoms in a crossed 
optical dipole trap (ODT) at 1064 nm, at a temperature 
that is twice the Rb condensation temperature Tc. The 
trap frequencies for Rb are 21 Hz in the horizontal (x, y) 
plane and 165 Hz in the vertical (z) direction; the trap 
frequencies for K are 1.37 times larger. The atoms are 
transferred into a 3D optical lattice in three steps. We 
first turn on a retro-reflected vertical beam in 150 ms to 
create a weak ID lattice. In the second step, the ODT is 
ramped off in 100 ms so that the two beams used for the 
ODT (which propagate along x and y) can be converted 
to lattice beams by allowing them to be retro- reflected. 
The intensities along all three directions are then ramped 
to their final values in 100 ms. The three lattice beams 
are derived from a common laser but individually fre- 
quency shifted. The x and y beams are elliptical with a 
200 X 40 jam waist and are linearly polarized orthogonal 
in the x-y plane; the z-beam has a circular waist of 250 
fim and is linearly polarized along x. We calibrate the 
lattice strength using Rb atoms with two different meth- 
ods (Kapitza-Dirac scattering pulse in a BEG [13] and 
parametric modulation of the lattice) and then account 
for the differences in mass and ac polarizability to deter- 
mine the lattice strength for KRb molecules. The values 
reported herein for the lattice depth in each direction are 
expressed in units of the molecule recoil energy, Er, and 
have an estimated 10% uncertainty. 

Once the atoms are loaded in the 3D lattice, we ramp 
an external magnetic field across an 5-wave Feshbach 
resonance at 546.78 G to form loosely bound "^^K^^Rb 
molecules with an efficiency of about 10%. With the 
Feshbach molecules at 5 = 545.8 G, where their bind- 
ing energy is /ix400 kHz, we use two-photon stimulated 
Raman adiabatic passage (STIRAP) to coherently trans- 
fer the Feshbach molecules to the ro-vibrational ground 
state [7], with a typical one-way transfer efficiency of 
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80%. All the molecules are in a single nuclear spin state 
in the rotational ground state, = O^ttin = O^mf = 
—A^mf^ = 1/2), following the notation defined in [8]. 
During this procedure, unpaired K and Rb atoms are 
removed using resonant light pulses. To measure the 
number of ground-state molecules in the lattice, we re- 
verse the STIRAP process and then image the resultant 
Feshbach molecules using absorption of a probe beam 
that is tuned to the imaging transition for K atoms. 
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FIG. 1: Loss of ground-state KRb molecules as a function 
of time in a 3D lattice with depths of 56, 56, and 70 Er in 
X, and respectively, where Er = h^k^/2m is the KRb 
recoil energy, k is the magnitude of the lattice beam wave 
vector, and m the molecular mass. Neglecting the very short 
time points (red solid circles), the number of molecules for 
times larger than 1 s (black solid circles) are fit to a single 
exponential decay, yielding a 1/e lifetime of 16.3±1.5 s. Inset: 
Lifetime in an isotropic lattice with a depth of 50 Er, with 
(blue open squares, 0.17 Debye) and without (black squares, 
Debye) an applied electric field. The lifetimes at 0.17 Debye 
(15di4 s) and Debye agree within uncertainty. 

Figure 1 shows a time-dependent evolution of the 
ground-state molecule population in the 3D lattice. In 
the first few lOO's of ms, the measured number of 
molecules exhibits relatively large variations in repeated 
iterations of the experiment and is consistent with some 
fast initial decay. In all our measurements of ground- 
state molecules in deep 3D lattices (for example, in the 
data for Fig. 2), we observe a similar feature. One 
possible explanation for this fast decay is collisions 
of the ground-state molecules with impurities, such as 
molecules in excited internal states that might be pro- 
duced in the STIRAP process. Fitting the data for times 
greater than 1 s to an exponential decay, which is consis- 
tent with a single-body loss mechanism, gives a 1/e life- 
time of 16.3±1.5 s. This is much longer than previously 
measured lifetimes of trapped ultracold polar molecules 
of about 1 s in an ODT [9] or in a ID lattice [11]. 

The long lifetime for ground-state molecules in a rea- 
sonably deep 3D lattice can be understood simply from 
the fact that the optical lattice localizes the molecules 
and therefore prevents bimolecular reactions. It was pre- 
viously seen that an applied electric field strongly in- 



creased the chemical reaction rate [10]. However, for 
molecules individually isolated in a 3D lattice, we ex- 
pect no dependence of the lifetime on the strength of an 
applied electric field. In the inset to Fig. 1, we show that 
indeed we do not observe any decrease of the lifetime for 
polarized molecules with an induced dipole moment of 
0.17 Debye. 

To understand what limits the lifetime of the molecules 
in the 3D lattice, we investigate its dependence on the 
lattice strength as summarized in Fig. 2. First, we ex- 
plore the transition from a 2D lattice (an array of one- 
dimensional tubes) to a 3D lattice. For a molecular gas 
confined in the tubes with no lattice in we find a life- 
time of ~1 s. However, as soon as a small lattice potential 
is added along the lifetime is dramatically increased, 
reaching 5 s at 12 Er and 20 s at 17Er (point a in Fig. 2). 
To verify that bimolecular reactions are the dominant 
loss mechanism, we have checked that the lifetime in un- 
corrugated tubes decreases significantly (to 0.1 s) when 
we apply an electric field (oriented along the tubes) that 
gives an induced dipole moment of 0.17 Debye. In addi- 
tion, the fact that we can place an upper limit of 10% of 
the initial number remaining at long times puts a limit 
on the contribution to our signal from tubes that are oc- 
cupied with only one molecule. 
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FIG. 2: Lifetime of KRb ground-state molecules in an optical 
lattice. Black circles: x and y lattice beams are fixed at 
56 Er per beam, while z is varied from to 136 Er (1 Er 
corresponds to a lattice intensity / = 0.025 kW/cm^). The 
lifetime reaches a maximum of 25d=5 s when the z lattice depth 
is 34 Er (point b). For higher lattice intensities, the lifetime 
decreases, which is consistent with loss due to off-resonant 
light scattering (dashed line). The open circles correspond 
to a 3D lattice where the radial confinement was also varied. 
The red squares correspond to lifetimes measured with an 
additional traveling- wave beam at 1064 nm illuminating the 
molecules in the 3D lattice. Point c (d) corresponds to the 3D 
lattice of point a with an intensity of 3.2 kW/cm^ (b with 3.7 
kW/cm^) plus the additional beam intensity of 2.3 kW/cm^ 
(3.5 kW/cm^). Solid lines: see text. 

We consider a number of factors (Supplementary In- 
formation) to understand the rapid suppression of loss as 
a function of the lattice strength in z. In general, Pauli 
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blocking for identical fermions and dissipation blockade 
effects (suppression of loss when the loss rate for particles 
on the same site is much larger than the tunneling rate 
[14]) can play a role in the lifetime of KRb molecules in an 
optical lattice. However, in the measurements reported 
here, the optical lattice is sparsely filled. For our identi- 
cal fermionic molecules, we expect Pauli blocking would 
give an even steeper function than we observe. Moreover, 
for a 5 Eji lattice in the lifetime of KRb molecules in 
the tube does not change significantly in the presence of 
an applied electric field. This observation suggests that 
an incoherent process, such as heating of the trapped 
gas, limits the lifetime in this regime of weak lattice con- 
finement. Instabilities in the optical phase of the lattice 
beams directly give rise to translational noise of the lat- 
tice, which can promote molecules to higher bands, where 
they have increased mobility and could then collide with 
other molecules. A simple theoretical model taking into 
account a constant heating rate, with collisions in higher 
bands happening on a timescale much shorter than the 
heating time, is consistent with the experimental obser- 
vation (red and blue solid lines in Fig. 2, with heating 
rates of 1 En/s (66 nK/s) and 2 En/s^ respectively). 

In Fig. 2, the lifetime reaches a maximum of 25±5 s 
indicated by point b. As the intensity is increased 
further, the lifetime starts to decrease, consistent with 
off-resonant photon scattering becoming the dominant 
loss mechanism. The rich internal state structure of 
molecules ensures that each off-resonant photon scatter- 
ing event has a high probability of causing the loss of a 
molecule from the ground state. To explore this effect, 
we added an additional traveling- wave beam with a wave- 
length of 1064 nm; this increases the photon scattering 
rate without increasing the trap depth and we observe 
a significant reduction of the lifetime due to the addi- 
tional light. We can extract the imaginary part of the 
polarizability of KRb molecules by fitting the lifetime as 
a function of the light intensity to l/{al). Here, a is the 
imaginary part of the polarizability at 1064 nm, which we 
determine to be (2.052 ± 0.009) x 10"^^ MHz/(W/cm2); 
this is consistent with a theory estimate for KRb [15, 16]. 

We have also explored the lifetime of KRb Feshbach 
molecules in the 3D optical lattice. It has been shown 
that these weakly bound molecules can be rapidly lost 
from an ODT due to collisions with atoms [17]. Even 
with removal of the Rb atoms and the RF transfer of 
the K atoms to a different hyperfine state, all previously 
measured lifetimes for KRb Feshbach molecules were less 
than 10 ms [18]. However, with the ability to create 
ground-state molecules, which do not scatter light that 
is resonant with the single-atom transitions, we can more 
efficiently use light pulses to remove any residual atoms. 
A few ms after the atom removal, we reverse the STI- 
RAP process and recreate a Feshbach molecule gas. For 
the case of molecules in the ODT (no lattice), we find 
that this extends the lifetime of the Feshbach molecules 
to 150 ms. Therefore, we can conclude that previous life- 
time measurements were likely limited by collisions with 



residual atoms. 

When we perform the procedure described above for 
KRb molecules in an optical lattice, we find that the puri- 
fied gas of Feshbach molecules can have a lifetime as long 
as 10 s. In Fig. 3 we explore the lifetime of the Feshbach 
molecules in the 3D lattice as a function of the magnetic- 
field detuning from the resonance {Bq = 546.78 G), where 
varying the magnetic field, changes the binding en- 
ergy and the size of the Feshbach molecules. We start by 
forming a purified sample of the Feshbach molecules in a 
strong 3D lattice with an intensity of 50 per beam at 
B = 545.8 G. The magnetic field is then ramped to its 
final value in 1 ms. At the end of the hold time in the 
3D lattice, B is ramped back to 545.8 G where we image 
the molecules. 
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FIG. 3: Lifetime of Feshbach molecules and confinement- 
induced molecules measured as a function of B. A purified 
sample of Feshbach molecules is held in an isotropic 3D opti- 
cal lattice (50 Er per beam, 20 kHz trap frequency). Near the 
Feshbach resonance, the loss rate due to photon scattering can 
be modeled (solid line) as a weighted sum of the free atom loss 
rate Tatom and a higher loss rate for tightly bound molecules 
-L molecule • The grey shaded area indicates the single atom 
lifetime, and its uncertainty, measured for the same experi- 
mental conditions. Inset: Lifetime of Feshbach molecules in a 
3D lattice as a function of the trap intensity, at 545.8 G (blue 
stars) and 543.18 G (green diamonds). The dashed and solid 
line fits are used to extract the scattering rates. 



Above the Feshbach resonance, the lattice potential al- 
lows for the existence of confinement- induced molecules 
that do not exist in free space [19]. We find that 
confinement-induced molecules have a lifetime (25 s) that 
is comparable to that for K or Rb atoms in the same 
trap. Below the Feshbach resonance, the molecule life- 
time decreases quickly when the magnetic field is ramped 
to lower values. Several Gauss below the resonance, the 
Feshbach molecule lifetime is reduced to ~1 s, which is 
still significantly longer than in the ODT. Overall, these 
results represent a two-orders-of-magnitude improvement 
in the lifetime of Feshbach molecules compared to a pre- 
vious measurement of KRb molecules in an optical lat- 
tice [20]. 

To understand the dependence of the lifetime on 5, 
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we can assume that the hfetime is hmited by off-resonant 
photon scattering from the lattice hght and consider two 
hmiting cases. For B ^ B^^ the photon scattering 
hmit is simply that for free atoms Vatom] for B <^ Bq, 
we have a higher photon scattering rate Vmoiecuie due 
to a larger wavefunction overlap with electronically ex- 
cited molecules [21]. In a two-channel model of the 
Feshbach resonance [22], the Feshbach molecule wave- 
function can be written as an amplitude Z^/^ times the 
bare "closed- channel" molecule wavefunction plus an am- 
plitude (1 — Z)^/^ times the "open channel" wavefunction 
that describes the scattering state of two free atoms. We 
then take the total photon scattering rate to be given 
by ZVmoiecuie + (1 " Z)Vatom' With pairs of atoms con- 
fined in an optical trap with a known depth, Z can be 
calculated straightforwardly with a coupled-channel the- 
ory [22, 23]. Using the measured loss rates for the limiting 
cases, Vatom and T molecule, this simple theory (sohd line 
in Fig. 3) without any additional adjustable parameters 
describes very well the experimental results (filled cir- 
cles). We note that the rate of atom- molecule collisions 
has been analyzed in a similar way [20, 24]. 

The assumption that the lifetime of the purified gas 
of Feshbach molecules in a 3D lattice is limited by only 
the photon scattering can be checked by varying the lat- 
tice beam intensity. This is shown in the inset of Fig. 3 
for two values of B. Similar to the case of ground- state 
molecules, we observe a rapid initial increase in the life- 
time going from no lattice (only the ODT) to a weak 
lattice. Following this initial rise, we observe a decrease 
in the Feshbach molecule lifetime as the lattice inten- 
sity is increased, consistent with loss due to off-resonant 
scattering of the lattice light. From the fits of Fig. 3 we 
extract an imaginary part of the Feshbach molecule's ac 
polarizability at 1064 nm of 15.9zbl.6 MHz/(W/cm2) for 
B = 545.8 G and 30±3 MHz/(W/cm2) for B = 543.18 G. 

In a simple model for the conversion of atoms to 
Feshbach molecules in a 3D lattice, one could assume 
100% conversion efficiency for individual lattice sites that 
are occupied by exactly one K atom and one Rb atom. 
Starting with a purified sample of Feshbach molecules 
prepared with round-trip STIRAP and atom removal as 
discussed above, we dissociate the molecules in the lat- 
tice by ramping B above Bq to 548.97 G. This should 
ideally produce only pre- formed pairs of atoms. We then 



ramp B back down to 545.8 G and measure the molecu- 
lar conversion efficiency. The result is (87 zb 13)%, where 
the uncertainty is dominated by fiuctuations in STIRAP 
efliciency in successive runs of the experiment. This high 
efficiency is far above the maximum of 25% observed in 
an ODT [25]; this indicates that optimizing the loading 
procedure in order to have a larger number of sites with 
exactly one Rb and one K atom would increase the overall 
conversion efficiency. For heteronuclear Bose-Fermi mix- 
tures, optimizing the number of pre- formed atom pairs in 
a lattice remains a challenge. Recent progress in this di- 
rection includes the characterization of a dual Bose-Fermi 
Mott insulator with two isotopes of Yb [26], as well as 
proposals to use interaction effects to optimize the lattice 
loading [27-29]. 

The capability demonstrated here for using a 3D lat- 
tice to freeze out chemical reactions, and thus prepare 
an ensemble of long-lived dipolar molecules, opens the 
door for studying many-body interactions in a gas of po- 
lar molecules in a lattice. For example, spectroscopy of 
rotational states of polar molecules in a lattice is a possi- 
ble approach to access correlation functions [30]. While 
individual molecules may experience long rotational co- 
herence times, dipolar interactions between neighboring 
lattice sites, which could have an interaction energy on 
the order of a few hundred Hz, will certainly require a 
systematic understanding of the many-body system to 
understand the resulting complex response. Future chal- 
lenges in studying these systems include achieving higher 
lattice filling factors, or correspondingly, lower entropy 
for a dipolar gas in a lattice. 
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I. LOSS RATES FROM GROUND-STATE MOLECULE COLLISIONS IN A CORRUGATED ID TUBE 

It was observed experiment ahy that the hfetime of molecules in an array of ID tubes is sensitive to the strength of 
an applied electric field. This electric field dependence is expected for collisions between identical dipolar fermions [1] , 
and hence implicates molecule- molecule collisions as the loss mechanism. It was also observed that, at zero electric 
field, the addition of a weak corrugating lattice along the tubes greatly enhances the molecule lifetime. In order to 
understand this effect, here we study lossy molecule-molecule collisions in a tube, at zero electric field and in the 
presence of a lattice along the z direction. 




FIG. SI: Loss processes within the lowest two bands of the z lattice, (a) Lowest band loss rate for two molecules on neighboring 
lattice sites (700). (b) Onsite loss rate for one molecule in the lowest band and one in the first excited-band on the same site 
(701). The offset A is in general present, due to the gaussian profile of the lattice beams. 

At low temperatures (/iK) and in zero electric field, a collision between two ground-state fermionic molecules 
in the same nuclear spin state occurs primarily in the p-wave channel. We model the lossy part of the molecule- 
molecule interaction with a p-wave pseudopotential, parametrized by a strictly imaginary scattering volume [2, 3]. 
The magnitude of the scattering volume, Vp ~ 350 nm^, is obtained by comparing ab initio calculations of the loss 
rate for two molecules in an isotropic 3D harmonic trap [4] to results using the pseudopotential. Assuming that most 
molecules are in the transverse ground-state of a tube, two representative loss processes for molecules in the same 
tube (on two adjacent lattice sites) are shown in Fig. SI. The loss rate for two nearest-neighbor molecules within 
the lowest band (700) is rapidly suppressed by Pauli blocking for increasing z-lattice depth. The loss rate for one 
molecule in the lowest band and one molecule in the first excited band on the same site (701) grows with increasing 
z-lattice depth, due to the increasing localization of the lowest band and first excited band Wannier orbitals {wo{z) 
and wi{z)^ respectively). 

We now address the question of whether the measured loss rates can be understood in terms of the fundamental 
two-body loss processes shown in Fig. SI, and their associated rates 700 and 701. Two molecules in the lowest band 
should be lost at a rate proportional to 700, where the constant of proportionality is determined in a semiclassical 
picture by the fraction of time those two molecules spend on neighboring sites. Since this fraction is clearly less than 
one — being driven down both by low density and by single particle localization due to the energy offsets between sites 
(A) — 700 is an upper bound to the actual loss rate of two molecules in the lowest band. As shown in Fig. S2, 700 is 
too small to explain the measured loss rates, so collisions between molecules within the lowest band are clearly not the 
loss mechanism. When one molecule is in the lowest band and one molecule is in the first excited band, the loss rate 
is proportional to 701, so long as hjoi ^ Ji, A is satisfied. It is still possible, however, that the loss rate can decrease 
with increasing lattice depths, since even as 701 is becoming larger the molecules are becoming more localized. When 
^7oi ^ Ji^^^ in the so-called Zeno suppression regime, the loss rate can actually be inversely proportional to 701. 
However, for the experimental data taken at low z-lattice depths (the first four data points in Fig. S2), the Zeno 
suppression scenario can be ruled out, and interband collisions should result in a loss rate that is proportional to 701. 



2 



7 ~ 7oi 



1000 



100 



10 



0.1 



0.01 



0.001 









Zeno suppression 




.»-»--"""■"""" 














O 










o 








o 




o 











■7 ~ </i/(7oi^^) 







10 20 30 

z-lattice depth (E' ) 



40 



FIG. S2: Relevant rates for molecules confined to the lowest two bands in the 2;-direction, as a function of the z-lattice depth. 
The black dashed-dotted curve is 700, the red solid curve is Ji/h, and the green dotted curve is 701. It is clear that 700 is too 
small to explain the measured loss rates (red circles, the error bars are smaller than the plot markers). The rate 701 increases 
with the 2;-lattice depth due to the increasing confinement of the lowest two Wannier orbit als, and eventually becomes larger 
than the first excited-band hopping rate (Ji/h) at a ^-lattice depth of around 23 ^r. For even larger 2;-lattice depths (green 
shaded region) one expects Zeno suppression to affect the loss rate of two molecules in different bands. 



This result is, however, inconsistent with the experimental observation that losses in corrugated tubes do not exhibit 
the electric field dependence of losses in uncorrugated tubes: When dipole-dipole interactions are taken into account, 
both 7oo and 701 should increase with electric field strength. The most likely explanation for this lack of electric field 
dependence is that, in the presence of the z-lattice, collisional physics does not determine the molecule lifetime, despite 
being ultimately responsible for the loss. Instead, a heating mechanism that promotes molecules from localized to 
mobile single-particle levels is what sets the time scale for loss. If the promotion rate into mobile states is slower than 
the collisional loss rate for mobile molecules, but faster than the collisional loss rate for localized molecules — which 
is suppressed by the diluteness of the molecular gas — then it will determine the molecule lifetime, thereby rendering 
the lifetime independent of electric field. 

To model this process phenomenologically, we assume that a heating mechanism imparts energy to the molecules 
at a fixed rate, and hence the promotion rate between successive bands is inversely proportional to the band gap 
between them. Microscopically, this picture is motivated by the following reasoning. Fluctuations in the transverse 
lattice beams — which provide a potential that is only approximately decoupled from the z-direction — should be the 
dominant source of heating at low z-lattice depths. While a detailed knowledge of these fluctuations is absent, one 
can in complete generality assume that to first order in the motion of a molecule along z couples to a time- 
dependent perturbation of the form X{z^t) = X±f{t)z^ where the intensity X± depends primarily on the transverse 
lattice strength. The resulting transition rate from a localized state \a) (Wannier orbital of the a^^ band along z) to 
a localized state |a + 1) (Wannier orbital of the next higher band along z) on the same lattice site is given by [5] 

r„,a+i = ^^^^^|f^Ka + lNa>|^ (SI) 

In the above, S{uJa) is the power spectrum of the noise /(t), evaluated at the transition frequency uj^ between \a) 
and \a + 1). For a power spectrum that is relatively flat over frequency ranges Auj ~ Er/H^ the matrix element 
|(a H- ^ gives the dominant scaling of the transition rate on the z-lattice depth, and thus we expect 

interband transitions of this variety to be suppressed by the band gaps. On a single site, if we take the energies of 
localized states to be the band centers eo^{v) (where v is the 2;-lattice depth measured in E^)^ then we can define 
interband transition rates Vcy^cyj^i{v) = £[ecx+i{v) — eoc{v)]~^ ^ where f is a constant heating rate. 

The loss rate for a given molecule is obtained by deciding on a localization criterion; for a molecule in the a*^ band 
sitting s sites from the trap center (in the z-direction), we consider that it is mobile whenever the bandwidth A Joe 



3 



provides enough kinetic energy for it to reach the trap center [6]. This criterion sets a critical lattice depth Vc{o^,s) 
in the a^^ band, below which a molecule at site s is mobile and above which it is localized. Our primary assumption 
is that the decay rate for a mobile molecule is much larger than the interband transition rate, in which case it can 
be chosen to be infinite. We then obtain the loss rate for a molecule on a given site and in a given band by solving 
coupled rate equations that describe the heating and loss. Because the experiments measure a molecule lifetime that 
is averaged over many molecules dispersed throughout the trap, we assume an initial distribution of molecules in the 
2:-direction that is gaussian with a radius of ^ = 14 lattice sites, which is based on experimental measurements 
of the molecule cloud, and average the results of our calculations over this distribution. The qualitative features 
of the results are not particularly sensitive to ^. We find that heating rates between £ = 1 E^/s (66 nK/s) and 
£ = 2 Er/s (132 nK/s) are consistent with the experimental data at low 2;-lattice depths, and molecule lifetimes as 
functions of the z-lattice depth for these two heating rates are shown in Fig. 2 of the manuscript. It should be noted 
that once the z lattice is sufficiently deep, and molecules remain localized for times comparable to the inelastic light 
scattering time, this model is no longer valid. 
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